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基于 地 物 高 光谱 和 无 人 机 多 光谱 的 黄河 三 角 洲 
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摘 要 : 土壤 盐 涡 化 是 限制 黄河 三 角 洲 地 区 农业 经 济 发 展 的 重要 因素 ,进一步 阻碍 了 农业 生产 。 为 了 探索 


无 人 机 影像 在 地 表 无 植被 覆盖 条 件 下 的 土壤 盐分 含量 


4 


反 演 状况 ， 以 黄河 三 角 洲 典型 区 域 为 研究 区 ， 获 取 地 


物 高 光谱 和 无 人 机 多 光谱 两 种 数据 源 与 样 点 土壤 盐分 含量 ， 通 过 优选 敏感 光谱 参量 ， 使 用 偏 最 小 二 乘 回归 
(Partial Least Squares Regression, PLSR) 和 随机 森林 (Random Forest, RF) 两 种 机 需 学 习 算 法 建立 土壤 盐 


分 含量 


间 的 敏感 性 最 高 ， 相 关系 数 为 -0.31。(2) 两 种 不 同 数 
(3) 基于 地 物 高 光谱 的 RF 模 型 (R?=0.54, RMSEv=3.30 g/kg) 优 于 基于 无 人 机 多 光谱 的 RE 模型 (R=0.54, 


KUE RMSRv=3.35 g/kg). (4) 结合 无 人 机 影像 采 


ES) 


演 模型 ， 实 现 研 究 区 的 土壤 盐分 含量 反 演 。 结 果 表 明 : (1) 高 光谱 1972 nm 波段 与 土壤 盐分 含量 


据 源 优化 后 的 RE 模型 均 优 于 PLSR， 且 稳定 性 更 好 。 


CH RF 模 型 对 研究 区 耕地 的 土壤 盐分 含量 进行 反 演 ， 


研究 区 总 体 以 轻 、 中 度 盐 渍 化 土壤 为 主 ， 对 作物 的 耕种 具有 一 定 程 度 的 限制 。 本 研究 构建 并 对 比 了 两 种 不 


同 源 数 据 的 黄河 三 角 洲 土壤 盐分 反 演 模型 ， 并 结合 


数据 源 的 优势 进行 优化 ， 探 索 了 地 表 无 植被 覆盖 情 


况 下 的 土壤 盐分 含量 反 演 方法 ， 对 更 精准 反 演 土地 盐 渍 化 程度 提供 了 参考 。 
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1 引 


土壤 盐 溃 化 是 世界 范围 内 重要 的 资源 及 生态 
问题 ,已 经 成 为 制约 区 域 农业 生产 和 经 济 可 持续 
发 展 的 主要 障碍 之 一 "”。 受 盐 汪 化 影响 的 土壤 
有 机 质 含量 较 少 、 肥 力 低 、 理 化 性 质 差 ， 作 物 不 


Il 
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A, 造成 了 土壤 中 次 生 盐 渍 化 严重 。 土 壤 退 化 、 
板结 ， 以 及 肥力 降低 等 问题 极 大 地 限制 了 农作物 
的 生长 ， 严 重 阻碍 了 黄河 三 角 洲 地 区 农业 经 济 的 
可 持续 性 发 展 ”。 

传统 的 土壤 盐 度 测量 方法 主要 通过 实地 采样 
和 实验 室 分 析 实 现 ， 费 时 费力 ， 很 难 对 农用 地 进 
行 大 面积 全 和 覆盖 的 动态 监测 “。 随 着 遥感 技术 的 
快速 发 展 ， 通 过 遥感 技术 对 土壤 盐分 含量 进行 准 
确 地 定量 反 演 成 为 了 许多 研究 者 关注 的 重点 “。 
在 多 源 影 像 方面 ，Wang 等 通过 构建 Cubist 模 
型 对 比 Landsat-8 OLI (Operational Land Imager) 
与 Sentinel-2 MSI (Multi-spectral Instrument) 传 
感 器 之 间 的 差异 ， 发 现 更 高 空间 分 辩 率 和 光谱 分 
辩 率 的 遥感 影像 有 助 于 提高 模型 反 演 的 准确 性 ; 
孙 亚 楠 等 ”通过 对 卫星 尺度 的 Landsat-8 OLI 和 
实测 高 光谱 数据 进行 融合 ， 构 建 内 蒙古 河套 灌区 
地 区 的 高 -多 光谱 融合 模型 ， 表 明 高 光谱 和 多 光 
谱 数 据 融合 能 够 提高 卫星 遥感 反 演 土 壤 盐分 含量 
的 精度 。Hu 等“ 使 用 电磁 感应 (Electromagnet- 
ic Induction, EMI) 设备 和 高 光谱 数据 构建 新 模 
型 ， 与 卫星 影像 构建 的 土壤 盐分 反 演 模型 作对 
比 ， 表 明 无 人 机 (Unmanned Aerial Vehicle, 
UAV) 多 光谱 数据 建 模 结果 更 优 。 搭 载 多 光谱 传 
感 器 的 无 人 机 在 农业 有 适 感 监测 方面 兼 具 多 方面 优 
势 ， 与 卫星 适 感 技术 相 比 ， 有 具有 高 空间 分 辨 率 、 
高 灵活 度 的 特点 ; 与 地 面 高 光谱 数据 相 比 ， 具 有 
图 谱 合 一 、 方 便 快捷 的 优势 。 在 地 表 植 被 长 势 与 
土壤 盐分 含量 间 的 联系 方面 ， 奚 雪 等 ”构建 无 
人 机 -卫星 一 体 化 的 偏 最 小 二 乘法 土壤 盐分 反 演 
模型 ， 在 黄河 三 角 洲 地 区 的 麦 区 构建 的 绿 度 指 数 
与 土壤 盐 渍 化 水 平 之 间 建 立 联系 ， 以 反 演 土壤 盐 
分 含量 。 贾 吉 超 等 '" 研究 了 冬小麦 种 植 区 与 土 


等 ' 引 入 红 边 波段 构建 新 的 光谱 指数 ,结果 表 
明 红 边 植 被 指数 模型 对 内 蒙古 灌区 的 土壤 盐分 含 
量 反 演 效果 更 优 。 多 光谱 影像 包含 地 物 的 光谱 信 
息 和 图 像 信息 ， 纹 理 特 征 是 图 像 空 间 灰 度 信息 的 
重要 属性 。 目 前 ， 纹 理 特 征 为 重要 的 影像 分 类 依 
Hk, TRE |) 利用 Sentinel-2 影像 结合 纹理 特征 
显著 提升 对 黄河 三 角 洲 地 区 中 度 盐 渍 土 的 分 类 效 
果 ， 有 效 提 高 影像 分 类 和 盐 渍 化 信息 监测 精度 。 
台 翔 ”通过 分 析 不 同 盐 分 条 件 下 土壤 和 植被 的 
光谱 特征 和 图 像 纹 理 ， 构 建 了 植被 覆盖 条 件 下 的 
土壤 盐分 反 演 模型 ， 加 入 纹理 特征 参数 能 有 效 提 
高 模型 精度 。 

目前 通过 遥感 技术 反 演 土壤 含 盐 量 的 研究 存 
在 以 下 两 个 问题 。 一 是 大 多 人 研究 集中 在 近 红 外 波 
段 附近 , 已 有 GF-5、HJ-1A 等 卫星 携带 了 高 光谱 
传感器 ， 但 受 限 于 空间 分 辩 率 和 波段 数目 ， 难 以 
在 较 小 尺度 内 精确 预测 土壤 盐分 含量 ; 二 是 随 着 
无 人 机 遥感 的 广泛 应 用 ， 纹 理 特 征 的 引入 能 和 否 提 
高 土壤 盐分 反 演 精度 。 因 此 ， 本 研究 通过 实地 采 
样 调研 ， 将 实测 的 地 物 高 光谱 与 无 人 机 多 光谱 的 
地 表 反 射 率 计 算 多 个 光谱 指数 并 进行 优选 ， 分 别 
建立 两 种 数据 源 的 土壤 盐分 反 演 模型 ， 并 结合 
自 数 据 特征 进行 优化 ， 对 比分 析 两 种 数据 源 的 反 
演 模 型 结果 并 评价 ， 最 终 使 用 无 人 机 影像 结合 精 
度 较 好 的 模型 反 演 研究 区 的 土壤 盐 渍 化 水 平 。 


2 数据 与 方法 


2.1 研究 区 概况 

人 研究 区 位 于 山东 省 东营 市 (118°07'— 
119°10'E, 36°55’—38°10'N), ， 是 黄河 三 角 洲 地 
区 典型 的 土壤 盐 渍 化 区 域 ， 属 于 暖 温带 半 湿 润 大 


培 盐 分 含量 间 的 联系 ， 发 现 冬小麦 种 植 范围 及 其 
变化 显著 受 土 壤 盐 分 状况 及 其 变化 的 影响 ， 土 壤 
含 盐 量 升 高 冬小麦 种 植 面积 又 减 。 黄 权 中 等 ” 
在 河套 灌区 进行 定点 监测 并 结合 OLI 影 像 ， 研 究 
发 现 土 壤 含 盐 量 与 玉米 、 向 日 获 的 长 势 和 产量 呈 
明显 负 相 关 。 在 传统 光谱 指数 优化 方面 ， 杨 宁 


陆 性 季风 气候 ， 年 均 降 水 量 少 但 蒸 散 量 较 大 ， 主 
要 土壤 类 型 为 潮 土 和 盐 土 ， 地 势 平坦 、 土 壤 盐 渍 
化 特征 明显 ”， 符 合 研 究 需要 。 综 合 东 营 市 区 域 
调查 结果 、 地 貌 土壤 以 及 农作物 的 分 布 情况 ,分 
别 在 研究 区 的 南部 和 东北 部 的 农作物 种 植 区 选择 
了 A、B、C 三 个 小 研究 区 ， 其 中 研究 区 A、B 主 
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要 种 植 作物 是 小 麦 ， 人 研究 区 C 分 为 距离 相近 C1 
和 C2 两 部 分 ，C1 种植 水 称 ，C2 为 含 盐 量 差异 性 
明显 的 裸 土 (图 1)。 研究 区 整体 分 布 上 具有 典型 


118°0'0"E 119°0'0"E 
a 1 


性 和 空间 异 质 性 ， 能 够 体现 黄河 三 角 洲 地 区 盐 渍 
化 土壤 的 分 布 趋势。 
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注 : 图 中 A、B、C 为 研究 区 


图 1 山东 省 东营 市 黄河 三 角 洲 试验 研究 区 分 布 图 


Fig. 1 Distribution map of Yellow River Delta experimental study area in Dongying city, Shandong province 


2.2 地 面 数据 采集 与 处 理 


于 2021 年 11 月 23 日 至 25 日 在 研究 区 内 进行 
了 实地 调查 ， 发 现 该 时 间 段 研究 区 降水 量 较 少 ， 
地 表 特 征 稳定 ， 土 壤 盐 涡 化 较为 明显 ， 便 于 提取 
土壤 端 元 的 光谱 特征 。 为 确保 均匀 随机 取 点 ,在 
研究 区 内 一 共 采 集 90 个 样本 点 ， 其 中 研究 区 A 
样本 10 个 ， 研 究 区 B 样 本 40 个 ， 研 究 区 C1 样本 
20 个 ， 研 究 区 C2 样 本 20 个 。 

使 用 华星 Al2 系列 全 球 导 航 卫 星系 统 
(Global Navigation Satellite System, GNSS) 接 
收 机 ， 测 量 每 个 样本 点 的 WGS84 大 地 坐标 。 待 
测 得 土壤 高 光谱 数据 后 ， 采 集 O~S cm 深度 的 土 
培 作为 待 测 盐 度 样本 。 每 个 样本 称 取 20 g 的 烘 干 
土壤 加 入 去 离子 水 ， 搅 拌 、 抽 滤 得 到 土壤 样本 
浸出 液 。 而 后 将 土壤 样本 浸出 液 倒 入 干燥 至 恒 
HANA ALP RETA, MATALA (30%) 
去 除 有 机 质 后 ， 放 入 烘箱 中 烘 干 至 恒 重 ,冷却 
后 称 量 ， 从 而 确定 每 个 样本 的 土壤 含 盐 量 (Soil 


Salinity, SS, g/kg)". 
2.3 光谱 数据 采集 及 预 处 理 


2.3.1 地 物 高 光谱 数据 采集 

在 研究 区 使 用 SVC HR-768i 型 便携 式 地 物 光 
41 (Spectra Vista Corporation， 纽 约 ， 美 国 ) 
测量 样本 点 表层 土壤 端 元 的 光谱 特征 。SVC HR- 
768i 型 地 物 光 谱 仪 共有 768 个 通道 用 以 接收 波长 
TE 350~2500 nm 之 间 的 电磁 辐射 数据 ， 在 350~ 
1000 nm 间 的 通道 带宽 小 于 1.5 nm, Æ 1000~ 
2500 nm 之 间 带 宽 为 4~8 nm。 观 测 时 天 气 晴 朗 ， 
传感器 探头 距离 地 面 1 m， 垂 直 正 射 样 点 土 体 自 
然 状态 的 表面 ， 每 次 测量 前 均 使 用 标准 白板 校准 
反射 率 ， 每 次 光谱 扫描 时 间 设 置 为 Ss， 同 时 对 
每 个 样 点 重复 测量 5 次 ， 取 反射 率 均值 作为 样 点 
的 土壤 高 光谱 反射 率 曲线 。 
2.3.2 无 人 机 多 光谱 数据 采集 

使 用 大 性 DJI 精 灵 4Pro 多 光谱 版 无 人 机 采集 
地 面 多 光谱 数据 ， 摄 像 头 带 有 1 个 可 见 光 成 像 的 
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彩色 传感器 和 5 个 多 光谱 成 像 的 单 色 传感器 : 蓝 
(B), & (G) 、 红 (R) 、 红 边 (RE) 、 近 红 
(NIR) 光谱 段 ， 中 心 波长 分 别 为 450、560、 
650、730 和 840 nm, ZG Tee a HL Bt F 
光谱 ， 近 红外 波段 通道 宽度 为 +26 nm， 其 他 波段 
均 为 +16 nm。 飞 行 前 ， 将 标准 灰 板 放置 在 研究 区 
内 对 传感器 进行 定 标 ， 试 验 期 间 无 人 机 飞行 时 间 
为 10:00 一 14:00， 飞 行 高 度 100 m， 影 像 采 集 间 
隔 2s， 晴 朗 无 风 。 

使 用 DJI Terra 软件 对 无 人 机 多 幅 影 像 进行 拼 
接 并 设置 地 理 参 考 ， 得 到 研究 区 的 各 波段 及 
RGB 影像 。 通 过 标准 灰 板 的 各 项 参数 对 多 波段 


影像 进行 辐射 定 标 ， 将 灰 度 级 转换 为 光谱 反射 
率 。 在 ENVI5.3 中 进行 无 人 机 影像 的 波段 融合 
和 裁剪 ， 结 合 实地 测量 样本 点 的 GPS (Global 
Positioning System) 坐标 数据 在 ArcGIS 10.6 中 提 
取样 本 坐标 点 对 应 像 元 的 5 个 波段 反射 率 数 据 。 


2.4 光谱 指数 的 选择 及 纹理 特征 的 引入 


2.4.1 光谱 指数 

依照 前 人 研究 进展 ， 计 算 多 个 光谱 指数 ， 如 
表 1 所 示 。 针 对 无 人 机 多 光谱 传感器 获取 的 红 边 
波段 影像 ， 引 入 使 用 红 边 波段 构建 的 新 型 光谱 
FR. 


RI 光谱 指数 计算 公式 


Table 1 Formula of spectral index 


光谱 指数 计算 公式 参考 文献 
归 一 化 植被 指数 (Normalized Difference Vegetation Index ,NDVI) NDVI = (NIR - R)/(NIR + R) (1) [16] 
土壤 调节 植被 指数 (Soil-Adjusted Vegetation Index , SAVI) SAVI = 1.5 X(NIR- R)/(NIR + R + 0.5) (2) [17] 
改进 非 线性 指数 (Modified Nonlinear Index, MNLI) MNLI = 1.5 X (NIR? - R)/(NIR? + R + 0.5) (3) [18] 
光谱 指数 SAVI*SR SAVI*SR = (NIR? - R)/[(NIR+R+0.5) XR] (4) [18] 
改进 叶绿素 吸收 比值 指数 (Modified Chlorophyll Absorption In MCARII =[(NIR - RE)- 0.2 X (NIR - R)](NIR/RE) 119] 
Reflectance Index] ,MCARI1 ) (5) 0 
改进 红 边 指数 (Modified Red Edge Simple Ratio, MRESR) MRESR = (NIR - R)/(RE - R) (6) [20] 
标准 化 红 边 指数 (Normalized Red Edge Index, NREI) NREI= RE/(NIR + RE + R) (7) [20] 
比值 植被 指数 (Ratio Vegetation Index , RVI) RVI = NIR/R (8) [21] 
最 优 植被 指数 (Optimal Vegetation Index, VI.,,,) VI = 1.45 X ((NIR? + 1)/(R + 0.45)) (9) [22] 
归 一 化 盐分 指数 (Normalized Differential Salinity Index ,NDSI) NDSI=(R- NIR)/(R + NIR) (10) [23] 
盐分 指数 SI1(Salinity Index) SH = SQRT(G X R) (11) [23] 
盐分 指数 SI3(Salinity Index) SI3 = SQRT(G2 + R°) (12) [23] 
盐分 指数 SIl-reg( Salinity Index) SIlreg = SQRT(G X RE) (13) [12] 
盐分 指数 SI3-reg( Salinity Index) SI3reg = SQRT (G° + RE?) (14) [12] 
soa T Based Normalized Differential NDSIreg = (RE - NIR)/(RE + NIR) (15) [12] 
FE: NIR、RE、R、G、B 分 别 代表 近 红 外 、 红 边 、 红 、 绿 、 蓝 波段 ， 在 高 光谱 数据 中 分 别 代 表 841、730、650、560 和 450 nm 波段 ; 


reg 代 表 计 算 中 用 红 边 波段 蔡 换 红 波 段 


2.4.2 ”纹理 特征 


关 性 ”。 由 于 灰 度 共生 矩阵 的 维度 较 大 ， 一 般 


影像 纹理 反映 了 土壤 表面 颜色 和 灰 度 的 变 
化 ， 这 些 变化 与 土壤 本 身 的 盐 度 、 含 水 量 等 因素 
BPA AMER SEF Sit IK EEE FEE 
法 提取 影像 纹理 特征 。 灰 度 共生 和 矩阵 (Gray Lev- 
el Co-Occurrence Matrix, GLCM) 是 一 种 在 特征 
提取 、 纹 理 分 析 以 及 影像 质量 评价 方面 广泛 应 用 
的 重要 方法 ， 用 以 描述 图 像 中 像素 灰 度 之 间 的 相 


不 直接 作为 区 分 纹理 的 特征 ， 而 是 基于 和 矩阵 构建 
的 统计 量 作 为 描述 图 像 纹 理 的 特征 。 

在 ENVI 5.3 中 对 原始 的 无 人 机 多 波段 影像 
进行 主 成 分 分 析 (Principal Component Analysis, 
PCA)， 对 多 波段 影像 进行 降 维 。 无 人 机 影像 的 
第 一 主 成 分 的 贡献 率 达 95% 以 上 ， 故 提取 影像 的 
第 一 主 成 分 ， 进 行 基于 灰 度 共生 和 矩阵 的 影像 纹理 
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特征 提取 ， 共 得 到 8 个 影像 的 纹理 特征 ， 包 括 均 
值 (Mean) 、 方 差 (Variance) 、 同 质 性 (Homo- 
geneity), 、 对 比 度 (Contrast) 、 差 异性 (Dissimi- 
larity) 、 (Entroy) 、 角 二 阶 窍 (Angular Sec- 


UAV 多 光谱 数据 


影像 拼接 、 辐 射 定 标 


CERD E 


(a) 土 壤 盐 分 含量 多 光谱 估 测 模型 


ond Moment, ASM) 和 相关 性 (Correlation) °”, 
25 土壤 盐分 含量 反 演 模型 的 构建 与 验证 


从 地 面 高 光谱 和 无 人 机 多 光谱 两 方面 分 别 构 
建 土壤 盐分 含量 模型 ， 思 路 如 图 2 所 示 。 


SVC 地 面 实测 高 光谱 数据 


反射 率 数 据 预 处 理 


光谱 参量 计算 j 波段 敏感 性 分 析 


(b) 土 壤 盐分 含量 高 光谱 估 测 模型 


图 2 土壤 盐分 含量 估 测 技术 流程 图 


Fig.2 Technical flow chart of soil salt content estimation 


(1) 样本 划分 。 将 90 个 样本 随机 分 成 60 个 
训练 样本 和 30 个 验证 样本 ， 对 60 个 训练 样本 的 
土壤 含 盐 量 数据 与 光谱 反射 率 数 据 计 算 的 多 个 光 
谱 指 数 进 行 Pearson 相关 性 分 析 。 

(2) 光谱 及 纹理 特征 筛选 。 无 人 机 多 光谱 数 
据 以 纹理 特征 和 光谱 指数 进行 筛选 ;地面 高 光谱 
在 原始 输入 波段 基础 上 ， 结 合 敏 感 波段 进行 进 一 
步 优化 。 根 据 土壤 的 光谱 特征 ， 使 用 无 人 机 采集 
的 蓝 (B)、 绿 (G)、 红 (R)、 红 边 (RE), if 
红外 (NIR) 共 5 个 波段 的 数据 计算 15 个 植被 指 
数 和 盐分 指数 ( 表 1)。SVC 地 物 高 光谱 仪 能 够 
采集 350~2500 nm 范围 内 多 个 窗 光 谱 段 的 反射 
率 数据 ， 但 是 与 无 人 机 多 光谱 相机 传感器 的 中 心 
波长 并 不 完全 重合 。 为 保持 一 致 性 ， 选 择 与 多 光 
谱 传 感 器 中 心 波长 尽 可 能 相近 的 光谱 段 。 选 择 实 
测 的 土壤 高 光谱 反射 率 数 据 中 450、560、650、 
730 和 841 nm 通道 分 别 代表 无 人 机 影像 中 的 B、 
G、R、RE、NIR 波段 。 特 征 优选 方法 采用 相关 
性 分 析 ， 排 除 无 法 通过 共 线 性 诊断 或 与 盐分 相关 
性 较 低 的 光谱 指数 ”。r 表 示 光 谱 参 量 间 的 相关 


系数 ,7 超过 0.9 时 认为 存在 共 线 性 问题 ， 在 0.8 
以 上 时 可 能 会 有 问题 O°, 

(3) 模型 选择 。 由 于 遥感 计算 中 各 光谱 指数 
之 间 存 在 比较 严重 的 多 重 相关 性 ， 为 减弱 模型 自 
变量 间 多 重 共 线性 问题 ， 采 用 偏 最 小 二 乘 回归 
(Partial Least Squares Regression, PLSR) 和 随机 
森林 (Random Forest, RF) 建 模 。 研 究 需 要 采 
用 多 个 光谱 指数 构建 模型 ， 偏 最 小 二 乘 回归 模型 
集合 了 主 成 分 分 析 、 线 性 回归 分 析 和 相关 分 析 的 
优势 ,在 预测 变量 间 存 在 高 度 相 关 甚 至 共 线 问题 
且 样 本 数量 较 小 时 ， 与 多 元 线性 回归 模型 
(Maximum Likelihood Estimation, MLR) 相 比 更 
具 优 势 。RF 模 型 将 多 个 互 不 联系 的 决策 树 构成 
森林 ， 对 多 元 共 线 性 问题 不 敏感 ， 结 果 对 样本 中 
的 缺失 数据 和 不 平衡 的 数据 比较 稳健 ， 是 机 器 学 
习 中 应 用 较为 广泛 的 算法 ”1。 选 择 与 土壤 盐分 
含量 显著 相关 的 光谱 参量 ,利用 Matlab R2021a 
中 构建 土壤 盐分 含量 的 估 测 模型 。 

(4) 精度 评价 。 以 土壤 含 盐 量 的 实测 值 与 预 
测 值 之 间 的 决定 系数 民 和 验证 集 均 方 根 误差 
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(Root Mean Squared Error of Validation, RMSEv) 
综合 评价 模型 的 训练 和 验证 精度 。 

(5) 盐 度 分 级 及 制图 。 选 择 较 优 的 一 种 建 模 
方法 建立 黄河 三 角 洲 土壤 盐分 含量 的 遥感 反 演 模 
型 ， 估 测 研 究 区 耕地 的 土壤 盐 渍 化 水 平 。 使 用 
ENVI 5.3 对 无 人 机 多 波段 影像 进行 处 理 ， 通 过 目 
视 解 译 提取 影像 中 耕作 区 域 。 综 合 对 比 所 构建 的 
PLSR 模 型 和 REF 模型 的 精度 ， 选 择 最 佳 的 裸 土 土 
培 盐 分 含量 遥感 估 测 模型 ， 用 其 对 研究 区 内 的 耕 
地 进行 无 人 机 遥感 尺度 的 土壤 盐分 含量 反 演 。 同 
时 对 耕地 的 土壤 盐 渍 化 程度 进行 分 级 ,将 土地 的 
盐 涡 化 程度 按照 国家 土壤 盐 渍 化 分 级 标准 划分 为 
5 个 等 级 ( 表 2)， 并 统计 面积 占 比 1。 

表 2 土壤 盐 渍 化 程度 分 级 标准 


Table 2 Grading standards for soil salinization 


th 4 RF Eh Eh Hjh 
土壤 盐 渍 化 等 级 Ea bss TL DR the 
mk me mikk mk 


土壤 含 盐 量 /(g' kg ) <1 1~2 2~4 4~6 >6 


3 ”模型 优化 与 结果 分 析 


3.1 土壤 盐分 与 光谱 特征 相关 性 分 析 


3.1.1 基于 地 物 高 光谱 的 土壤 盐分 指数 分 析 及 
优化 

对 比 低 盐 度 与 高 盐 度 样本 光谱 特征 差异 ， 分 
析 地 面 多 光谱 数据 各 波段 反射 率 与 土壤 盐分 含量 
之 间 相 关 性 (图 3) 结果 表明 ， 对 350~2500 nm 
间 全 部 的 光谱 段 进 行 Pearson 相关 性 分 析 ， 在 
350~1350 nm 间 为 正 相 关 ，1350~2500 nm 间 为 
HIK; 在 1400 和 1850 nm 波段 附近 存在 剧烈 
ZSAE, 1972~1978 nm 间 的 反射 率 对 盐 度 的 敏感 
性 最 高 ， 其 中 1972 nm 波段 的 相关 系数 ”为 最 敏 
感 的 -0.31。 可 见 光 波段 均 未 达 极 显著 相关 ， 绿 
光 和 蓝光 相关 性 最 强 ， 红 边 波段 敏感 性 较 弱 r 接 
WE 0。 短 波 红 外 (Short Wavelength Infra-Red, 
SWIR) 与 近 红 外 波段 相 较 ， 具 有 更 加 敏感 的 负 
相关 性 ， 考 虑 使 用 SWIR 构建 光谱 参量 。 

分 析 地 面 高 光谱 数据 与 土壤 含 盐 量 相 关 性 


加 
一 一 高 盐 度 样本 光谱 反射 率 


局 


光谱 反射 率 
L 


1 ! -0.4 
1000 2000 3000 
波长 /nm 
图 3 各 光谱 段 与 土壤 含 盐 量 间 的 相关 性 分 析 
Fig. 3 Correlation analysis between spectral bands and soil 


salt content 


(图 4 (a) )， 结 果 表 明 ， 基 于 近 红外 波段 构建 的 
36 16 48 Be, MCARII, NREI, MRESR, NDSI- 
reg 达 到 0.05 KÆ E m 5 4H; NREI 和 MRESR 
两 个 光谱 指数 与 土壤 盐分 含量 达到 显著 相关 ， 且 
彼此 间 不 存在 极 强 的 共 线 性 问题 MRESR 为 正 
相关 ，NREI 为 负 相 关 。 通 过 全 光谱 段 的 相关 性 
分 析 ，R 和 NIR 波 段 的 反射 率 敏 感性 较 差 ， 选 择 
与 土壤 盐分 含量 相关 性 更 优 的 波段 重新 计算 光 
谱 参 量 以 提高 模型 精度 ， 所 以 选择 使 用 SWIR 
1972 nm 代替 原始 光谱 指数 中 的 NIR 841 nm。 结 
果 表 明 ， 使 用 SWIR 1972 nm 重 构 光 谱 参 量 后 ， 
所 有 包含 近 红 外 波段 运算 的 光谱 指数 与 土壤 含 盐 
量 间 的 相关 性 均 得 到 了 大 幅度 提升 (图 4 (b) )。 
除 VI, 及 SI 类 指数 外 ， 其 他 光谱 参量 均 达 极 显著 
相关 。 在 相关 系数 的 矩阵 中 可 以 看 出 ， 光 谱 参 量 
间 的 共 线 性 问题 有 所 加 重 ，MCARI]1 与 MRESR 
之 间 的 7 为 0.91， 认 定 两 个 指数 间 存 在 共 线 性 问 
题 ， 因 此 ， 后 续 建 模 中 选择 MNLI、SAVI*SR、 
MRESR 和 NDSI 作 为 构建 模型 的 自 变量 。 
3.1.2 基于 无 人 机 多 光谱 的 土壤 盐分 相关 性 
分 析 

无 人 机 影像 提取 地 面 点 反射 率 数 据 ， 与 土壤 
盐分 相关 性 分 析 ， 筛 选 满 足 0.01 水 平 极 显 著 相 关 
且 不 存在 严重 共 线 性 问题 的 光谱 指数 (图 5)。 
SAVI*SR、VI SIL, SI3, SIl-reg 和 SI3-reg 与 
样本 含 盐 量 均 达 到 极 显 著 相 关 。 在 盐分 指数 中 ， 
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(a) 基 于 近 红 外 波段 构建 的 高 光谱 指数 
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(b) 基 于 短波 红外 构建 的 高 光谱 指数 


注 :** 表 示 在 0.01 水 平 ( 双 侧 ) 上 显著 相关 ;* 表 示 在 0.05 水 平 ( 双 侧 ) 上 显著 相关 
图 4 高 光谱 指数 与 土壤 盐分 相关 性 分 析 


Fig. 4 Correlation analysis between hyperspectral indices and soil salinity 


BR NDSIŻh, SIL, SI3, SIl-reg, SI3-reg 敏感 性 
均 表现 较 好 ， 但 是 指数 间 存 在 极 强 的 共 线 性 问 
题 ， 相 关系 数 + 接 近 1。 相 较 SWIR， 可 见 光 区 间 
内 敏感 性 较 低 ， 红 边 波段 不 是 土壤 盐 度 的 敏感 波 
段 区 间 ， 绿 光波 段 的 贡献 性 更 高 。VIL, 与 4 个 盐 
分 指数 均 纯 在 较 高 的 相关 性 ， 存 在 共 线 性 问题 ， 
选择 VIF<10 的 Vl (r= -0.35), SAVI*SR (r= 
0.35), SI3-reg (r= 0.42) 三 个 光谱 指数 构建 无 
人 机 尺度 的 土壤 盐分 估 测 模型 。 


Ss Ss vs. Šo iS eS 
& SF LF LELE SSS o KIS 

SS | s 
NDVI 
SAVI 
MNLI 
SAVI*SR 
MCARII 
MRESR 
NREI 
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Vlopt 
NDSI 
si 
SI3 
SI1-reg 
SI3-reg 
NDSI-reg 


注 : ** 表 示 在 0.01 水 平 ( 双 侧 ) 上 显著 相关 ;* 表 示 在 0.05 水 平 ( 双 
侧 ) 上 显著 相关 
图 5 无 人 机 多 光谱 指数 与 土壤 盐分 相关 性 分 析 
Fig. 5 Correlation analysis between multi-spectral indices 


of UAV and soil salinity 


3.2 无 人 机 影像 的 纹理 特征 分 析 


考虑 无 人 机 多 光谱 影像 数据 的 优势 ， 引 入 影 
像 的 纹理 特征 作为 新 的 模型 自 变 量 。 对 PCA 处 
理 后 的 影像 计算 ， 获 得 8 个 纹理 特征 波段 后 ， 按 
照 坐 标 提取 各 样本 点 对 应 的 纹理 特征 并 与 盐 度 进 
行 敏感 性 分 析 ( 表 3)。 无 人 机 影像 提取 纹理 特征 
会 造成 的 数据 元 余 ， 对 8 个 纹理 特征 两 两 之 间 计 
算 相关 系数 r+-， 排 除 呈 高 度 相关 的 特征 变量 , 减 
弱 数 据 宛 余 对 模型 精度 的 影响 。 同 质 性 (x = 
0.16) 、 对 比 度 (r= -0.15) 和 差异 性 (x = 
-0.16) 3 个 特征 向 量 的 敏感 性 接近 ， 在 无 人 机 图 
像 中 各 样本 点 的 区 分 度 不 高 ， 对 盐 度 的 敏感 性 有 
限 。 由 于 纹理 特征 受到 图 像 噪声 的 一 定 影响 ， 在 
SPAHR TEP, tHe MAM C= 
-0.26) 、 角 二 阶 矩 (r = 0.27) 两 个 特征 向 量 满 
足 在 0.05 水 平 上 显著 相关 ， 其 他 特征 向 量 的 相关 
性 均 表 现 不 高 ， 对 模型 精度 的 提高 没有 意义 。 但 
E, MAMA ZME AAFEERI, r 
为 0.97， 认 定 二 者 间 存 在 严重 共 线 问题 ， 存 在 严 
重 共 线性 问题 的 纹理 特征 对 于 土壤 盐 度 反 演 模型 
的 构建 是 没有 帮助 的 ， 因 此 本 研究 选择 角 二 阶 矩 
与 其 他 光谱 参量 作为 构建 无 人 机 多 光谱 反 演 模型 


的 自 变量 。 
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68 
表 3 无 人 机 影像 各 纹理 特征 与 土壤 盐分 含量 相关 性 
Table 3 Correlation between texture characteristics and soil salt content 
SS Mean Variance Homogeneity Contrast Dissimilarity Entroy ASM Correlation 
Ss 1 
Mean 0.19 1 
Variance -0.18 0.26" 1 
Homogeneity 0.16 -0.44" -0.60" 1 
Contrast -0.15 0.33” 0.83” -0.81" 1 
Dissimilarity -0.16 0.40" 0.74” -0.96" 0.94” 1 
Entroy -0.26" 0.23" 0.52" -0.74" 0.52" 0.67" 1 
ASM 0.277 -0.19 -0.43" 0.69" -0.43" -0.59" -0.97" 1 
Correlation  -0.06 -0.09 0.19 0.04 0.06 0.01 0.14 -0.09 1 


注 : ** 表示 在 0.01 水 平 ( 双 侧 ) 上 显著 相关 ; * 表 示 在 0.05 水 平 ( 双 侧 ) 上 显著 相关 


3.3 基于 地 面 高 光谱 数据 的 土壤 盐分 估 测 
模型 的 优化 与 验证 


使 用 高 光谱 波段 优化 后 的 4 个 光谱 指数 作为 
模型 的 自 变量 ,分 别 构建 PLSR 和 RF 土壤 盐分 
估 测 模型 (图 6) 。 通 过 优化 前 后 的 模型 精度 
对 比 ，PLSR 模型 的 RMSEv 从 3.59 g/kg 减少 到 
3.38 g/kg，R* 从 0.11 提 升 到 0.20; RF 模型 的 RM- 
SEv 从 4.30 g/kg 减少 到 3.30 g/kg， 提 升 更 加 明 
显 ， 有 较 好 的 精度 ， 尼 从 0.37 提 升 到 0.54。 可 以 
A tt} SWIR 对 裸 土 土壤 盐分 含量 反 演 精度 的 提升 


25 训练 集 n=60 


是 有 效 的 ， 土 壤 中 盐分 含量 的 差异 在 短波 红外 波 
段 的 相关 性 变化 更 具 特 征 性 。 从 图 6 可 以 看 出 ， 
改进 后 的 两 种 模型 对 裸 土 土壤 盐分 含量 的 反 演 能 
力 依旧 有 限 ，RF 模 型 的 精度 要 优 于 PLSR 模 型 ， 
RF 模型 对 土壤 盐分 反 演 的 效果 更 优 ， 但 是 受到 
土壤 样本 盐分 含量 分 布 不 均匀 的 限制 ， 对 盐分 含 
量 过 高 的 土壤 样本 很 难 做 到 准确 地 反 演 ， 并 且 对 
一 些 低 盐 度 土 壤 的 反 演 存在 一 定 高 估 ， 但 是 对 大 
部 分 中 、 低 盐 度 样本 拟 合 趋势 明显 ， 对 高 盐 度 样 
本 的 预测 结果 也 有 一 定 程 度 改善 ， 具 有 更 好 的 拟 
合 精度 。 


验证 集 n=30 
一 一 实测 值 
—— RF-SWIR 
一 一 PLSR-SWIR 


样本 数 /个 
图 6 波段 优化 后 偏 最 小 二 乘 回归 模型 (PLSR-SWIR ) 和 随机 森林 模型 (RF-SWIR ) 土 壤 盐 分 估 测 结果 
Fig. 6 Soil salinity estimation results of band optimized spectral indices for PLSR model(PLSR-SWIR )and 
RF model(RF-SWIR ) 
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3.4 基于 无 人 机 多 光谱 影像 的 土壤 盐分 反 
演 模型 的 构建 和 验证 


将 角 二 阶 矩 作为 自 变 量 带 入 土壤 盐分 反 演 模 
型 ， 与 无 纹理 特征 的 模型 相 比 〈 表 4) ， 模 型 的 尼 
和 RMSEv 都 有 所 提升 。 结 合 无 人 机 纹理 特征 
PLSR 模 型 的 RR 为 0.26，RMSEv 为 3.59 g/kg, RF 
模型 RR 为 0.54，RMSEv 为 3.35 g/kg, HAR 
纹理 特征 有 助 于 提升 无 人 机 遥感 反 演 土壤 盐分 含 
量 的 精度 。 相 较 地 面 高 光谱 数据 模型 (RMSEv= 
3.30 g/kg) ， 高 光谱 模型 更 为 稳定 ， 拟 合 的 效果 


25 训练 集 n=60 


| 


含 盐 量 /(g-kg 1!) 


更 好 。 同 样 ， 多 光谱 RF 模型 的 估 测 效果 优 于 
PLSR 模型 ， 对 佑 测 高 盐 度 的 样本 准确 性 有 了 一 
定 提 升 ， 可 以 看 出 对 验证 集 预测 效果 有 了 一 定 提 
Fr (图 7)。 受 到 土壤 盐 度 样本 和 算法 模型 自身 限 
制 ， 不 能 对 超越 样本 范围 的 数据 进行 预测 ， 也 难 
以 对 盐分 含量 过 高 和 过 低 的 土壤 样本 准确 反 演 。 
对 高 盐 度 土壤 的 反 演 准确 度 低 ， 是 造成 模型 整体 
精度 较 差 的 主要 原因 ， 但 是 盐 度 大 于 15 g/kg 样 
本 的 反 演 结果 符合 盐 土 的 评价 标准 ， 在 无 人 机 影 
像 制图 中 影响 较 小 。 


验证 集 n=30 
一 一 实测 值 
—=—RF-T 
] —+— PLSR-T 


图 7 结合 纹理 特征 与 无 人 机 多 光谱 指数 的 偏 最 小 二 乘 回归 模型 (PLSR-T) 和 随机 森林 模型 (RF-T) 土 壤 盐 分 估 测 结果 
Fig. 7 Soil salinity estimation results of PLSR(PLSR-T )and RF(RF-T)models combining texture features and 


UAV based spectral indices 


3.5 基于 无 人 机 多 光谱 影像 的 研究 区 土壤 
盐分 含量 反 演 


通过 PLSR 和 REF 模 型 的 对 比 ， 选 择 使 用 无 
人 机 多 光谱 数据 构建 佑 测 土壤 盐分 的 RF 模型 。 
按照 土壤 盐 汗 化 分 级 标准 将 研究 区 内 的 土壤 盐 
分 含量 划分 为 5 个 等 级 ( 表 2) ， 统 计 不 同 盐 渍 
化 等 级 土壤 的 面积 所 占 总 面积 的 百分比 ， 评 价 耕 
地 质量 。 从 研究 区 内 各 等 级 盐 渍 化 土壤 的 面积 统 
计 情 况 〈 表 5$) 可 以 看 出 的 研究 区 C 与 A、B 相 
较 ， 整 体 的 土壤 盐 渍 化 问题 更 加 严重 ， 土 壤 盐 渍 
化 程度 最 高 ， 与 实地 调查 情况 相符 ， 重 度 盐 渍 化 
和 盐 土 的 总 面积 超过 了 30%。 和 研究 区 A、B 均 为 


轻 、 中 度 盐 渍 化 耕地 ， 研 究 区 C 为 中 、 重 度 盐 渍 
化 。 结 合 实际 地 理 空间 的 分 布 ， 符 合 越 靠 近 黄河 
入 海口 ， 土 壤 盐 渍 化 问题 更 严峻 的 事实 。 由 于 土 
坏 盐 渍 化 严重 影响 作物 生长 ， 研 究 区 A、B 耕 种 
的 主要 作物 是 小 麦 ， 而 研究 区 C 耕 作 的 是 水 稳 。 


4 讨论 


在 基于 遥感 技术 对 区 域 土壤 盐分 含量 定量 反 
演 的 研究 中 ， 对 地 表 裸 土 盐分 含量 的 准确 分 析 相 
对 匮乏 ， 仍 然 需要 深入 齐 析 和 探索 。 本 研究 使 用 
SVC 地 物 光 谱 测量 仪 获取 地 面 高 光谱 数据 ， 对 样 
本 点 进行 多 次 测量 取 均值 ， 保 证 土壤 高 光谱 数据 
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表 4 土壤 盐分 含量 反 演 模型 对 比 引入 短波 红外 波段 后 的 光谱 指数 与 土壤 盐分 
Table 4 Inversion model of soil salinity content 含量 的 相关 性 在 波长 1900 nm 附近 显著 提升 ， 这 
精度 评价 与 前 Bx LS 中 盐 碱 士 光 谱 反 射 率 在 该 波段 阶 
bra Boa : 与 前 人 研究 FP RC EET De LE 
R RMSEv(gke) 近 的 陡 降 相 吻 合 。 在 土壤 盐 涡 化 过 程 中 产生 Na- 
Jonu J Cl 等 化 合 物 ， 这 些 化 学 物质 与 水 在 键 态 结构 上 的 
PLSR-SWIR 0.20 3.38 xŠ i Fd 
svc HIE Li Ay ap E SE A IRB MER EO? HE 
RF-NIR 0.37 4.30 ‘ 
1400~1500 nm, 1800~1900 nm 形成 敏感 性 较 强 
RF-SWIR 0.54 3.30 rr er ern) ee 
= ae 的 光谱 豚 ， 对 土 六 盐分 含量 间 的 差异 分 析 具 有 更 
WR A PR ab 
ME 90 as 多 深入 探索 的 可 能 。 
av 2 7 a ee 本 研究 通过 构建 高 光谱 和 多 光谱 两 种 数据 
RF-T 0.54 3.35 类 型 的 裸 土 土壤 盐分 含量 估 测 模型 ， 并 进行 优 


注 : PLSR-NIR, RF-NIR 为 近 红 外 高 光谱 模型 ，PLSR-SWIR、 
RF-SWIR 为 短波 红外 高 光谱 模型 ; PLSR、RF 为 未 引入 纹理 特 
征 的 多 光谱 模型 ，PLSR-T、RF-T 为 结合 纹理 特征 的 多 光谱 模型 


A B Cl 
0 30m 0 100 200m 
Lt | 
nace | 4 
Soil salinity grade [| | 2 | | 5 
(g/kg) 0 100 m 
E ; — 


图 8 研究 区 土壤 盐分 反 演 图 
Fig. 8 Inversion map of soil salinity in the study area 
表 5 黄河 三 角 洲 研究 区 土壤 盐 渍 化 等 级 面积 占 比 统计 
Table 5 Proportion statistics of soil salinization grade area 


in the study area 


pA 


人 


单位 : % 
E 等 级 
研究 区 

1 2 3 4 5 

0 43.1 54.2 1.5 1.2 
B 0 25.9 70.2 2.6 13 
Cl 0 10.4 58.7 19.4 11.5 
C2 0 4.3 26.1 28.6 41.0 


特征 的 准确 提取 。 使 用 无 人 机 获取 多 光谱 数据 ， 
对 研究 区 先进 行 无 人 机 测量 ， 后 进行 实地 高 光谱 
测量 ， 保 证 无 人 机 影像 的 纹理 特征 不 受 人 为 因素 
的 干扰 ， 尽 可 能 保持 土壤 的 自然 状态 ， 保 证 反 演 
结果 更 加 准确 。 


化 和 分 析 。 优 化 后 ， 两 种 数据 源 的 PLSR 模型 
分 别提 升 验 证 集 RMSEv 分 别 优化 到 3.38 g/kg 
和 3.59 g/kg，RF 模 型 的 验证 RMSEv 分 别 优 化 到 
3.30 g/kg 和 3.35 g/kg， 表 明 无 人 机 影像 结合 图 谱 
言 息 能 够 充分 发 挥 自 身 优势 。RF 在 泛 化 能 力 和 
平衡 误差 方面 表现 良好 ， 但 受 限于 机 器 学 习 的 鲁 
棒 性 ， 存 在 一 定 的 过 拟 合 问题 ETA 
光谱 数据 的 REF 反 滨 模型 (R=0.54， 验 证 RMSEv 
二 3.30 g/kg) 与 基于 无 人 机 多 光谱 的 RE 模型 (R? 
二 0.54， 验 证 RMSEv 二 3.35 g/kg) 相 较 ， 精 度 较 
高 。 曹 肖 奕 等 优选 地 面 高 光谱 指数 反 演 土壤 
电导 率 ， 研 究 表明 ， 土 壤 电导 率 的 敏感 波段 位 于 
2000 nm 附近; 机 器 学 习 算 法 能 够 显著 提高 土壤 
电导 率 的 估算 精度 。 本 研究 结果 与 该 结论 一 致 。 

纹理 特征 是 反 演 土壤 盐分 含量 的 重要 依据 。 
研究 选用 描述 影像 灰 度 分 布 一 致 性 与 匀 质 性 的 角 
二 阶 矩 对 土壤 盐 度 进行 区 分 呈 。 灰 度 分 布 均匀 
时 ， 角 二 阶 矩 较 大 ; 灰 度 均匀 性 弱 时 ， 角 二 阶 和 矩 
较 小 。 因 影像 为 无 植被 覆盖 的 裸 土 ， 表 面 均匀 性 
强 的 盐 土 和 均匀 性 弱 的 非 盐 碱土 存在 明显 差异 ， 
角 二 阶 矩 较 好 区 分 这 种 差异 。 无 人 机 多 光谱 影像 
的 纹理 特征 受到 多 方面 的 影响 ， 与 点 对 点 采集 的 
地 面 多 光谱 数据 相 比 ， 纹 理 特 征 不 能 直观 表现 出 
含 盐 量 对 土壤 的 影响 程度 ， 而 地 面 高 光谱 能 够 直 
接 表现 出 盐分 差异 引起 的 光谱 特征 变化 ， 使 得 无 
人 机 模型 的 验证 精度 低 于 地 面 高 光谱 模型 有待 
开展 对 影像 纹理 特征 的 深入 探究 和 应 用 。 盐 碱土 
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板结 程度 不 同 使 无 人 机 影像 纹理 产生 差异 ， 对 土 
培 盐 渍 化 程度 的 判断 具有 一 定 提升 。 土 壤 的 光谱 
反射 率 不 仅 与 盐分 含量 相关 ， 土 壤 光 谱 特 征 受 含 
水 量 、 有 机 质 和 颗粒 度 等 复杂 理化 性 质 的 影 
响 5， 在 裸 土 土 壤 盐 度 的 反 演 方面 仍 需要 进 一 
步 研 究 。 土 壤 盐 度 在 部 分 地 块 差异 性 较 强 、 变 化 
性 大 ， 例 如 研究 区 C， 对 于 盐 度 在 15 g/kg 以 上 的 
样本 反 演 结果 偏 低 ， 对 模型 整体 精度 的 影响 
较 大 。 


5 结 论 


本 研究 基于 地 面 高 光谱 数据 和 无 人 机 多 光谱 
数据 提出 了 构建 裸 土 土壤 盐 度 的 估 测 模型 。 选 择 
统计 分 析 法 和 机 器 学 习 法 分 别 构建 PLSR 和 REF 的 
裸 土 土壤 盐分 含量 反 演 模型 ， 并 对 模型 进行 验证 
和 优选 。 获 得 以 下 主要 结论 。 

(1) REF 盐分 反 演 模型 均 优 于 PLSR 模 型 ， 其 
中 最 优 的 是 高 光谱 RF 土壤 盐分 含量 反 演 模型 (R= 
0.54, RMSEv=3.30g/kg) 。 短 波 红外 优化 光谱 参 
量 有 助 于 提升 估 测 土壤 盐分 含量 的 精度 。 

(2) 无 人 机 影像 纹理 特征 的 引入 ， 图 谱 信 息 
的 结合 进一步 提升 土壤 盐分 含量 反 演 精度 ， 但 无 
人 机 模型 验证 精度 低 于 高 光谱 模型 。 

(3) 黄河 三 角 洲 地 区 土壤 盐 渍 化 水 平 在 地 理 
区 位 上 差异 明显 ， 南 部 地 区 的 耕地 盐 汗 化 水 平 较 
低 ， 适 宜 种 植 耐 盐 碱 弱 的 作物 ;靠近 东北 部 黄河 
入 海口 土壤 盐 涡 化 水 平 较 高 ， 宜 种 植 耐 盐 碱 较 强 
的 作物 。 
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Machine Learning Inversion Model of Soil Salinity in the 
Yellow River Delta Based on Field Hyperspectral and 
UAV Multispectral Data 
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Abstract: Soil salinization in the Yellow River Delta is a difficult and miscellaneous disease to restrict the development of agri- 
cultural economy, and further hinders agricultural production. To explore the retrieval of soil salt content from remote sensing 
images under the condition of no vegetation coverage, the typical area of the Yellow River Delta was taken as the study area to 
obtain the hyperspectral of surface features, the multispectral of UAVs and the soil salt content of sample points. Three represen- 
tative experimental areas with flat terrain and obvious soil salinization characteristics were set up in the study area, and 90 sam- 
ples were collected in total. By optimizing the sensitive spectral parameters, machine learning algorithms of partial least squares 
regression (PLSR) and random forest (RF) for inversion of soil salt content were used in the study area. The results showed that: 
(1) Hyperspectral band of 1972 nm had the highest sensitivity to soil salt content, with correlation r of -0.31. The optimized 
spectral parameters of shortwave infrared can improve the accuracy of estimating soil salt content. (2) RF model optimized by 
two different data sources had better stability than PLSR model. RF model performed well in terms of generalization ability and 
balance error, but it had some over-fitting problems. (3) RF model based on ground feature hyperspectral (R’*=0.54, verified 
RMSE=3.30 g/kg) was superior to the random forest model based on UAV multispectral (R’=0.54, verified RMSE=3.35 g/kg). 
The combination of image texture features improved the estimation accuracy of multispectral model, but the verification accura- 
cy was still lower than that of hyperspectral model. (4) Soil salt content based on UAV multi-spectral imageries and RF model 
was mapped in the study area. This study demonstrates that the level of soil salinization in the Yellow River Delta region is sig- 
nificantly different in geographical location. The cultivated land in the study area is mainly light and moderate salinized soil 
with has certain restrictions on crop cultivation. Areas with low soil salt content are suitable for planting crops in low salinity 
fields, and farmland with high soil salt content is suitable for planting crops with high salinity tolerance. This study constructed 
and compared the soil salinity inversion models of the Yellow River Delta from two different sources of data, optimized them 
based on the advantages of each data source, explored the inversion of soil salinity content without vegetation coverage, and can 
provide a reference for more accurate inversion of land salinization. 

Key words: soil salinity; remote sensing; ground feature hyperspectral; UAV multispectral; partial least squares regression; ran- 


dom forest; machine learning 
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